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Abstract 

The balanced vehicular traffic model is a macroscopic model for vehicular traffic 
flow. We use this model to study the traffic dynamics at highway bottlenecks either 
caused by the restriction of the number of lanes or by on-ramps or off-ramps. The 
coupling conditions for the Riemann problem of the system are applied in order 
to treat the interface between different road sections consistently. Our numerical 
simulations show the appearance of synchronized flow at highway bottlenecks. 
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Highway bottleneck 
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1 Introduction 



The balanced vehicular traffic model (BVT model), which was first intro- 
duced in [10], generalizes the macroscopic traffic model of Aw, Rascle [1,13] 
and Greenberg [3] by introducing an effective relaxation coefficient into the 
momentum equation of traffic flow. This effective relaxation coefficient can 
become negative, resulting in multivalued fundamental diagrams in the con- 
gested regime. Such negative effective relaxation coefficients follow from fi- 



Email addresses: {f .siebel, w.mauser}@iggf .geo.uni-muenchen.de (Florian 
Siebel and Wolfram Mauser), {salissou, rascle}@math. unice.fr (Salissou 
Moutari and Michel Rascle). 



Preprint submitted to Elsevier Science 



2 February 2008 



nite reaction and relaxation times of drivers (see [10]). For related ideas, see 
Greenberg-Klar-Rascle [5], Greenberg [4]. 

In a previous work [11] we studied the behavior of the BVT model at a bot- 
tleneck. There, we manipulated the partial differential equations describing 
traffic flow at the bottleneck by artificially resetting the average velocity in 
order to model a speed restriction. The results obtained in that study showed 
the basic behavior of the BVT model and its potential to explain the observed 
patterns of traffic flow [9]. Although the procedure of resetting quantities, i.e. 
the average velocity in that case, was locally restricted, it leaves the question 
whether the model can adequately describe synchronized flow at a bottleneck. 
In order to show this, the current paper systematically studies the traffic dy- 
namics at highway bottlenecks in the BVT model by numerical means. To 
do this in agreement with the underlying partial differential equations we use 
the coupling conditions for the Riemann problems at the interface between 
different highway sections. We focus the discussion on two setups: In the first 
setup we study a bottleneck caused by the narrowing of a highway from three 
lanes to two lanes. In the second setup we study a two-lane highway with an 
on-ramp and an off-ramp. 

The outline of the paper is as follows. In Section 2 we summarize the theory 
of the coupling conditions for the Riemann problem at intersections. Whereas 
Section 3 presents numerical results for a highway bottleneck caused by the 
restriction of the number of lanes from three lanes to two lanes, Section 4 
presents the numerical results for a two-lane highway with an on-ramp and an 
off-ramp. In Section 5 we introduce the possible generalization to an arbitrary 
junction. We finally summarize our results in Section 6. 



2 Coupling conditions 



Before we have a closer look at the coupling conditions, let us repeat the prin- 
cipal equations of the BVT model. The evolution equations for the density p 
of vehicles and the average velocity v are described by the following hyperbolic 
system of balance laws 



a W v-u M) + 3{ P v {v -u {p))) =bMp(u(p) _ vy (2) 

Here, u(p) denotes the equilibrium velocity, therefore the expression v — u(p) = 
w describes a distance to equilibrium. The quantity b(p, v) is the effective relax- 
ation coefficient. For the case where b(p,v) becomes negative, there are addi- 
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tional equilibrium velocities, i.e. high-flow branch and the jam line. Whereas 
the high-flow branch is metastable for intermediate densities and unstable 
for high densities, the jam line is unstable for intermediate densities and 
metastable for high densities. For a detailed discussion see [11]. Note that 
the pseudo-momentum p(y — u(p)) is not conserved due to the non- vanishing 
term on the right-hand side of Eqn. (2). This source term plays an essential 
role for the traffic dynamics on road sections, but it is neglected for the situ- 
ation where one is interested in the Riemann problems at intersections, since 
it is never a delta-function. 



2. 1 Background 

Piccoli and Garavello [2] appear to have been the first to propose an inter- 
section modeling by using the Aw-Rascle "second order" model of traffic flow 
[1]. In their approach, only the mass flux is conserved but not the pseudo- 
momentum. In [8], Herty and Rascle proposed another approach in which 
mass flux and pseudo-momentum are both conserved. But they maximized 
the mass fluxes at the intersection with some arbitrary given homogenization 
coefficients. In [7], the latter approach has been generalized by maximizing 
the total mass flux at the junction without fixing any condition. In fact, the 
homogenization coefficients are not arbitrary but obtained directly from the 
mass flux maximization. Another approach also based on the mass flux max- 
imization is given in [6]. Our approach in the current paper is similar to the 
latter one, in particular we maximize the total flux at the junction, at the 
same time conserving the pseudo-momentum of the original "Aw-Rascle" sys- 
tem. Here, we are particularly concerned with the BVT model [10]. We note 
that our treatment of the homogenization problem which naturally arises in 
a merge junction is different from the one in [8], [7], see Remark 1. 




a) b) c) 

Fig. 1. Definition of the quantities at junctions, which are necessary to construct the 
boundary fluxes. The junction in a) consists of one incoming and one outgoing road. 
The panels b) and c) display a merge junction and a diverge junction, respectively. 
The direction of the flow is from left to right. 
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2.2 Junction consisting of one incoming and one outgoing road 



In the following we restrict the discussion to the Riemann problem at the in- 
terface between two road sections. The road section 1 is located upstream and 
the road section 2 downstream of the interface, where the flux functions have 
to be determined. Let £/f = (pi,PiVi) be the state vector at the interface 
upstream in the road section 1 and let = (pt , pt v t) be the state vector at 
the interface downstream in the road section 2 (see panel a) of Fig. 1). For each 
road section j we introduce the function Wj of the state variable U = (p, pv) 

Wj(U) =v -Uj{p). (3) 

According to Herty and Rascle [8] the fluxes at the interface between the two 
road sections (see Eqs. (l)-(2)) 



f = q 



i N 



(4) 



Mur); 

can be calculated from the expression 

g = min(di(pf),s 2 (/4))> ( 5 ) 

where the demand and supply functions d\{p) and s 2 (p) and the density p\ 
are defined below. 

The density p\ is defined as the intersection point of the two curves v — v%, 
W2(U) = Wi(Ui) and can be obtained by solving the implicit equation 

u 2 (pl)-v£ + Wl (Ur) = 0. (6) 

Let us define the function 

Vdi(p) = pui(p) + pwi(E/T). (7) 

For a monotonously decreasing, differentiable function U\(p) the function %i(p) 
has a single maximum at location p 1; which can be obtained by solving the 
implicit equation 

^(pO = pVi(pi) + + w!(yr) = o. (8) 
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It is now possible to define the demand function 1 

4W = H?: (f) 

[%i(Pi), if P> Pi- 
Let us define the function 

Vs2{p) = pMp) + pwi{Ui). (10) 

For a monotonously decreasing, differentiable function u 2 (p) the function f] S 2{p) 
has a single maximum at location p 2 , which is determined by solving the im- 
plicit equation 

Vaifc) = P2u' 2 (p 2 ) + u 2 (p 2 ) + wi(C/f) = 0. (11) 
With this function we define the supply function as 

Mp) = ("*ff (12) 
[77 s2 (p), if p > p2- 

We numerically solve the implicit Eqs. (8) and (11) using the method of nested 
intervals. With the expression for the fluxes at the interface between the two 
road sections (see Eqn. (4)), we obtain the necessary boundary values at the 
interface for the conservative update scheme described in detail in [10]. 

2. 3 Merge junction 



We are now interested in the boundary fluxes at a merge junction, see panel 
b) of Fig. I. In order to define the demand functions on the incoming road 
sections i — 1, 2 we first define the functions r)di(p) as follows, 

Vdi(p) = pUi{p) + pWi{U~). (13) 

Let us denote the maxima of the corresponding curves as pV We can then 
define the demand function of the road section i as 

[VdiiPi), if P > Pi- 

In order to define the supply function in the road section 3, we first introduce 
quantities fa describing the fraction of cars entering from the road section % 



1 Note that in [12] we exchanged the terms demand and supply in comparison to 
the standard notation used here. The demand describes the maximum flow that the 
road section 1 can deliver to the road section 2. 
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into the road section 3. When we assume that the incoming fluxes passing 
through the junction are proportional to the incoming demands, we have 



ft = ^%=-y (15) 

With this definition, we follow [6] and do not consider these fractions as part 
of the optimization problem as in [7]. We define the homogenized value w 3 for 
the quantity w defined in Eqn. (3), i.e. 



w* 3 ^J2P jWj (U-), (16) 
j'=i 

and with the latter the function 

Vss(p) = pu 3 (p) + pw* 3 (17) 
which reaches its maximum value at p 3 . 

Remark 1 Here, we assume that the velocity on road 3, near the junction, 
is given by v = u 3 (p) + w 3 . This is in contrast with [8] and [7], where the 
mixture of cars from both incoming roads 1 and 2 is assumed to produce an 
homogenized flow, with a nonlinear relation between p and v which expresses 
that the cars from roads 1 and 2 microscopically share the available space. The 
resulting mixture rule is more appropriately described in Lagrangian (mass) 
coordinates, but is definitely not in the above form. 

In other words, the assumptions in [8] and [7] are incompatible with the (sim- 
pler) formula (17), which we assume here. Naturally, in practice the difference 
is not necessarily significant. We will come back to this point elsewhere. 

We are now able to define the supply function 

s ,M = Hf « p l h ' (18) 

[rj s3 (p) : if p > p 3 - 

This function will be evaluated below at the intersection point of the two 
curves v — v 3 and w 3 {U) = w 3 , which corresponds to a density p\ fulfilling 
the implicit equation 

uM)-vt + w* 3 = 0. (19) 
Finally, we define the downstream boundary fluxes in the road section i as 



f~ = q(3i 



( 1 \ 

\wi(Ur) J 



(20) 
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and the upstream boundary fluxes in the road section 3 as 



\ w 3 



where 



q = min (J2 d j(Pj )> s 3(p 



(21) 



(22) 



Note that the above boundary fluxes are conserved through the intersection, 
and are bounded from above by the demands and supplies. 



2.4 Diverge junction 



In this section we study the boundary fluxes at a diverge junction, see panel 
c) of Fig. 1. Again we set 

Vds{p) = pu 3 {p) + pw 3 (U^) (23) 

and denote the maximum of that function as p 3 . The demand function in the 
road section 3 is defined as 

<h(p) = H p }: (24) 

[Vd3{p3), A P > P3- 

For the definition of the supply function we first have to prescribe the fractions 
of cars intending to enter from road section 3 into road section 1 and road 
section 2, which we describe with the quantities a 3 \ and 0:32. We stress that 
these quantities - in contrast to [6] - need not agree with the actual percentages 
of the flow from road section 3 entering into road section 1 and road section 
2, see Eqs. (31)-(34) below. Under the assumption that all cars remain in the 
network, i.e. a 3 i + a 32 = 1, we can set 



"31 = «, (25) 

«32 = (1"«), (26) 

with a G [0, 1]. For the outgoing roads sections k — 1, 2 we define 

Vsk = pu k (p) + pw 3 (Uz) (27) 

and denote the maxima of these functions The supply function on the 
outgoing road section k then reads 

/ v jvsk(pk), if p < Pfc, , , 

s k{p) = < . . _ (28) 
\Vsk{p), np>pk- 
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We also have to determine the densities p\, at which these supply functions are 
evaluated. These densities are calculated from the intersection of the curves 
v — v£ and Wk(U) = w 3 (U 3 ), which reduces to the implicit equations 

u k (pl)-4 + ^(U^)=0. (29) 

Finally, we define downstream boundary fluxes of road section 3 as 

• ( i ^ 

f 3 — q 

and the upstream boundary fluxes in road section 1 as 

/ 1 

K ™z(U 3 ) 



and in the road section 2 as 



f-2 



\ 
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where 



Mu 3 ) 

q 1 = min (ad 3 (p^ ), Si(pJ)) , 
q 2 = min ^(1 - a)d 3 (p 3 ), s 2 (p5)) » 



and 



(30) 



(31) 



(32) 



(33) 



(34) 



g = ?i + g 2 - (35) 
Note again that the above boundary fluxes are conserved through the interface, 
and are bounded from above by the demands and supplies. 



3 Lane reduction on a highway 



We study the traffic dynamics for the setup depicted in Fig. 2. The highway 
under study consists of two 7 km long road sections. The road section 1 con- 
sists of three lanes whereas the road section 2 consists of two lanes. Note that 
in the mathematical description the transition from two to three lane is im- 
mediate, the length of the merging segments is neglected. As in [11] we use 
the equilibrium velocity function of Newell 
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Fig. 2. Sketch of the highway under study. The highway consists of two road sections 
of 7 km length. The road section 1 consists of three lanes whereas the road section 
2 consists of two lanes. We use periodic boundary conditions, i.e. the road section 
1 is also located downstream of the road section 2. 

with parameter values u m = 160km/h, A = 3600 [1/h/lane], p m = 160 [1 /km/lane] 
and an effective relaxation coefficient 



b(p,v) 



u—v 



if (3(p,v)(u(p) - v)—a o >0, 
ifP(p,v)(u{p)-v)-d c <0, 
P(p,v), else, 



d c 
u—v 



(37) 



TUr. 



u{p) — v + aiAt>| + a 2 Av 



and 



A-u(p) = tanh 



o 3 — )(u(p) 

Pm ' 



CPm 



1 1 

P Pm 
0.1 S, CLi - 



(3* 



(39) 



with parameters a c = 2 m/s 2 , d c = —5 m/s 2 , T = 0.1 s, ai = —0.2, a 2 = —0.8, 
03 = 7 and c = — 14 km/h. Thus the maximum density of the road section 1 is 
480 vehicles/km. The road section 2 can support 320 vehicles/km at maximum. 
The initial data to start the numerical simulations consists of equilibrium 
data on the two road sections. We prescribe a constant vehicle density po in 
both road sections, setting the initial velocity to v = u(p ). We choose the 
constant po to be independent of the number of the lanes, the corresponding 
scaled densities in each road section follow from dividing po by the number 
of lanes of that road section. In the following we perform a parameter study 
of po, varying the quantity between 50 [1/km] and 300 [1/km] in steps of 50 
[1 /km] . Figure 3 displays the simulation results for the density (left column) 
and velocity (right column) for simulations covering two hours. Note that, 
although the initial data are in equilibrium in each road section, the coupling 
conditions at the interface between the two road sections do not guarantee the 
equilibrium during the evolution. 

For a density po = 50 [1/km], a small region of higher density and lower 
velocity forms between about 5.5 km and 7 km. This region corresponds to 
data located in the fundamental diagram on and scattered around the jam 
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line. Clearly, this congested region is fixed at the bottleneck and therefore 
cannot correspond to a wide moving jam. The structure is supported by the 
bottleneck, i.e. the insufficient capacity of the road section 2 to carry the 
corresponding free flow rates in the road section 1. This region corresponds to 
synchronized flow. For a density po = 100 [1/km], the dynamics becomes more 
complicated, but finally a synchronized flow region of extended width (ranging 
from about 2 km to 7 km) forms. Only in a small region of the road section 1 
(between km and 2 km) traffic is in free flow. For a density p — 150 [1/km], 
the synchronized flow region covers the entire road section 1. Moreover, a 
wide moving jam travels through the two road sections, see the bright pink 
structure in the velocity plot. Increasing the density still further leads to even 
wider wide moving jams. The velocities inside these jams decrease with the 
increase of the initial density p . For p — 300 [1/km] velocities of less than 1 
[km/h] are reached inside the wide moving jam. Note that the wide moving 
jams are not affected by the interfaces between the two road sections of the 
highway. They travel upstream with an almost constant speed. 



4 Bottlenecks caused by on-ramps and off-ramps 



In the second setup we analyze by numerical means a two-lane highway with an 
on-ramp and an off-ramp. The simulation setup is displayed in Fig. 4. For our 
simulations we chose a length of 7 km for two-lane road sections 1 and 3 each, 
and a length of 10 km for the one-lane road section 2. For the parameterization 
of the equilibrium velocity curve and the effective relaxation coefficient, we 
use again the values given in Eqs. (36)-(39). We start the simulations with a 
constant density in equilibrium on all road sections of po — 50 [1/km/lane] 
and vary the percentage of cars a aiming to enter the road section 1 from the 
road section 3, see Eqn. (25). The numerical results are summarized in Fig. 5. 
For small to intermediate values of the parameter a, the simulations develop 
a state, where the road section 1 is almost empty, i.e. traffic is in free flow in 
the road section 1. This can be easily understood by realizing that for small 
values of a most of the cars use the by-pass road section 2. In contrast, traffic 
is in the congested regime in the road section 3. For sufficiently high values 
of a (see the results for a = 0.5) synchronized flow develops in front of the 
on-ramp in the road section 1. For a = 0.7 there is only a small region of free 
flow remaining in road section 1 which disappears after a time of about 0.8 h. 
For a parameter value a = 0.9 a region of narrow moving jams forms in front 
of the off-ramp in the road section 3. 
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5 Extension to a general junction 



For a given junction n, let us denote by 5~ and <5+, respectively the set of all 
incoming roads to n (indexed % in the sequel) and the set of all the outgoing 
roads to n (indexed k in the following). We require the equations (l)-(2) to 
hold on each road of 5~ U <5+. The percentage of cars on the road i intending 
to go to the road k are denoted by an-, such that Vi G 5~, J2k a ik = 1- These 
coefficients are assumed to be known. 

Let U~ = (Pi,p~v~), V« G 5~ and Uj~ = (pt 7 Pk v k)> V A; G <5+, respectively 
the boundary values on the incoming and outgoing roads. We denote by (3 ik , 
such that VA; G S^,J2iPik = 1, the proportion of cars on the road k coming 
from the road i. We set 



Put = = aikMp p ( , yieS-,yke5t (40) 



with the demand functions di as defined in Eqs. (13)-(14). The homogenized 
w on the outgoing roads near the junction are as follows 

< = E PikMUr), VA; G fi+. (41) 

With these quantities we define the supply functions as in Eqs. (17)-(18) 
for arbitrary k G 5+. For all k G 5+, the intermediate state of density p k 
on the outgoing road is given by the intersection point between the curves 
v k (U) = v£ and w k (U) = w* k . 

To obtain the flux q on each road one has to solve the following maximization 
problem 



max ^2 Qk subject to 


(42a) 


kes+ 




0<qi< di(p-),Vi G 5~; 


(42b) 


0<q k <s k (pl),Vke5+; 


(42c) 


Qi = E Pikq k ,Vi G 5~; 


(42d) 


kes+ 




Qk< J2 aikdi(p~),Vk G £+. 


(42e) 
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6 Conclusion 

We have studied the balanced vehicular traffic model at highway bottlenecks 
caused by the reduction of the number of lanes and the effects of an on-ramp 
and an off-ramp. To this aim we performed numerical simulations changing 
the initial density and the routing parameter at the off-ramp respectively. For 
the lane reduction setup the numerical results show that already for moderate 
densities, the synchronized flow forms at the bottleneck. The width of the 
synchronized flow region increases with increasing density. For large densities, 
wide moving jams appear which travel with a constant velocity upstream. 
Wide moving jams are not affected by the interface between the highway 
sections, i.e. by a change in the number of lanes on the highway. For the 
setup with an on-ramp and an off-ramp our numerical simulation show that 
synchronized flow can form upstream of the on-ramp, but also in front of an 
off-ramp, where narrow moving jams can emerge. 

The theory of the coupling conditions described in this paper can be applied 
to the balanced vehicular traffic model at a general junction, thus guarantee- 
ing the conservation of the fluxes in the corresponding Riemann problems at 
intersections. 
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Fig. 3. Traffic dynamics at the bottleneck caused by the reduction of the number 
of lanes on a highway. The column on the left shows the evolution of the vehicle 
density in units [1/km/lane], the column on the right the corresponding evolution of 
the velocity in units [km/h]. The different rows correspond to different simulation 
runs varying the initial density p = po as indicated. See the text for a detailed 
description. 
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Fig. 4. Sketch of the simulation setup. We study the dynamics on a two-lane highway 
with an on-ramp between the road sections 1 and 3 and an off-ramp between the 
road sections 3 and 1. On- and off-ramp form the boundaries of the one-lane road 
section 2. 
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Fig. 5. Traffic dynamics at highway bottlenecks caused by on-ramps and off-ramps. 
The plot shows the traffic dynamics on the two-lane highway, the road section 1 
corresponds to the region between and 7 km, the road section 3 corresponds to 
the region between 7 and 14 km. We vary the percentage of cars intending to enter 
from the road section 3 into the road section 1, i.e. the value of a, in the range 
between 0.1 and 0.9 (different rows). The column on the left shows the evolution of 
the vehicle density in units [1/km/lane], the column on the right the corresponding 
evolution of the velocity in units [km/h]. 
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